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The Friedberg-Lee model is studied at finite temperature and density. By using the finite temper- 
ature field theory, the effective potential of the Friedberg-Lee model and the bag constant B{T) and 
B{T, jj.) have been calculated at different temperatures and densities. It is shown that there is a crit- 
ical temperature Tc — 106.6MeV when jj, — OMeV and a critical chemical potential /j. ~ 223.1MeV 
for fixing the temperature at T = 50MeV. We also calculate the soliton solutions of the Friedberg- 
Lee model at finite temperature and density. It turns out that when T < Tc (or fj, < nc), there 
is a bag constant B{T) (or B{T,fi)) and the soliton solutions are stable. However, when T > Tc 
(or /i > /ic) the bag constant B{T) = OMeV (or B{T,jj,) — OMeV) and there is no soliton solution 
anymore, therefore, the confinement of quarks disappears quickly. 
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I. INTRODUCTION 



^ . . . 

O i' It is widely believed that at sufficiently high temperatures and densities there is a deconfinenient phase transition 

JL , between normal nuclear matter and quark-gluon plasma (QGP), where quarks and gluons are no longer confined in 

ij ' hadrons[l] [2] . The theoretical and experimental investigation of QGP is one of the most challenging problems in high 

r^ I energy physics. The main goal of heavy ion experiments is to create such form of matter and study its properties. 

Due to the property of confinement and asymptotic freedom, QCD can not yet be used to study low-energy nuclear 

physics. The challenge to nuclear physicists is to find models which can bridge the gap between the fundamental 

theory and our wealth of knowledge about low energy phenomenology. Some of these models have been proved to be 

(^ ' successful in reproducing different properties of hadrons, nuclear matter and quark matter [3, 4, 5, 6, 7, 8], among 

T^ ' which Friedberg-Lee model[9], also referred to as nontopological soliton model[10], has been widely discussed in the 

^^ ■ past decades, see for instance Refs.[ll][12][13]. The model consists of quark fields interacting with a phenomenological 

'^ [ scalar field a. The a field is introduced to describe the complicated nonperturbative features of the QCD vacuum. It 

shows an intuitive mechanism for the deconfinenient phase transition. In physical vacuum state, the physical value of 

a is large and quark mass is more than IGeV which makes it energetically unfavorable for the quark to exist freely, 

C ■ . so that the effective heavy quarks have to be confined in hadron bags. But with the temperature and/or density 

^^ ' of the system increasing, the physical value of a decreases and in turn the effective quark mass drops down, the 

thermodynamic motion leads to a deconfinement of the effective light quarks. 

The properties of Friedberg-Lee model, mostly in the mean field approximation, are widely studied by Wilets and 

. , his co-workers[10][12]. The model has been very successful in describing the static properties of isolated hadrons. 

S . Moreover, during the past several years, Friedberg-Lee model was extended to finite temperature and density to 

study deconfinement phase transitions[14][15][16][17][18]. In the framework of finite temperature field theory, most 

of those studies focus on the effective potential of the Friedberg-Lee model at finite temperature and density, but 

seldom considered the value of the soliton configurations at high temperature and density. In Our previous work[19] 

the improved quark mass density-dependent model(IQMDD) was studied at finite temperature. By using the finite 

temperature field theory, the effective potential of the IQMDD model and the bag constant B{T) have been calculated 

at different temperatures, and it is shown that there is a critical temperature at T ~ llOMel^ when the deconfinement 

phase transitions take place. In the meantime, the soliton solutions at different temperatures are obtained. The 

Lagrangian of the Friedberg-Lee model we use is almost identical to that of IQMDD model which was studied in 

our previous work, especially their potentials take the same forms, but there are two basic differences between the 

Friedberg-Lee model and IQMDD model suggested first in Ref.[20]. First, for the the Friedberg-Lee model, inside the 

bag (a region of the metastable vacuum), the quarks have zero rest mass and only kinetic energy, on the contrary, 

according to the IQMDD model, the masses of u, d quarks are given by niq = -^—{q = u, d, u,d), B is the vacuum 
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energy density inside the Bag, and ns is the baryon number density. Second, in the Friedberg-Lee model, the a field is 
a phenomenological scalar field which is introduced to describe the complicated nonperturbative features of the QCD 
vacuum, but for the IQMDD model, the <t field is the meson field, and the coupling of quarks and the meson field a 
is introduced to improve the QMDD model [8] by including the quark-quark interactions. So that the parameters of 
the Friedberg-Lce model are significantly different from the ones in the IQMDD model, but the discussions adopted 
in the previous one can still be applied. 

In the resent years, the research of the the deconfinement phase transition at finite temperature and density has 
received a great deal of attention, then it is interesting to extend the Friedberg-Lec model to the finite temperature 
and density to study the detailed behaviors of the soliton solutions at different temperatures and densities. Since 
there are lots of works on the effective potential at finite temperature, in present paper, we especially focus on the 
study of the value of the soliton configurations at different temperatures and densities, in doing so we can give more 
detailed information of the the deconfinement phase transition. 

The organization of this paper is as follows. In the following section we review the Friedberg-Lee model and our 
selected parameters. In Sec. Ill, we give detailed calculation of the effective potential of the Friedberg-Lee model and 
the bag constant as functions of temperature and density. The soliton solutions of the Friedberg-Lee model at different 
temperatures and densities are presented in Sec. IV, while in the last section wc present our summary and discussion. 

II. THE FRIEDBERG-LEE MODEL 

The Lagrangian of the Friedberg-Lce model has the form[9] 



C = ij{id-gc7)iP + -d^ad^a-Uia), (1) 



which describes the interaction of the spin-^ quark fields ip and the phenomenological scalar field a with the coupling 
constant g, U{a) is a potential describing the nonlinear self-interactions of the a field. We will discuss only u and d 
quark in this paper. The potential for the a field is chosen as 

t/(-) = 1-^ + ^-^ + 5-^ + 5, (2) 

b^ > iac (3) 

The model parameters a, b and c are fixed so that U{ij) has a local minimum at (Tq ~ 0, and a absolute minimum at 
a large value of the a field 



2c 



8ac 



(4) 



The condition (3) ensures that the absolute minimum of U{a) \s aX a — a^, ^ Q. ctq represents a metastable vacuum 
where the condensates vanishes, a^ corresponds to the physical or nonperturbative vacuum. The difference in the 
potential of the two vacuum states is the bag constant B. If we take U{a^) = 0, the bag constant B can be expressed 



as 



The Euler-Lagrange euqations corresponding to (1) are given by 

iid-ga)ij = 0, (6) 

aa+-^+glp^p = 0. (7) 

da 

In the mean-field approximation (MFA), the scalar field a is taken as a time- independent classical c-number field, and 
we only consider a fixed occupation number of valence quarks (3 quarks for nucleons, and a quark-antiquark pair for 
mesons). Quantum fiuctuation of the bosons and effects of the quark Dirac sea are thus to be neglected. 

In the following, by solving the Euler-Lagrange euqations, the numerical calculation indicates that the a field has 
a soliton solution which is of a spherical cavity-like structure, the metastable vacuum is inside the cavity while the 



physical vacuum is outside, and the quarks are confined inside the cavity, however, with increase of temperature or 
density, the solutions of the equation of motions will become a damping oscillation, i.e. the soliton solutions are 
melted away and the deconfinement phase transitions take place. 

In the spherical case, the a field is spherically symmetric, and valence quarks are in the lowest s-wave level. Then 
the scalar field a and the Dirac equation functions can be written as 

(T(r,t) = a(r), (8) 

^{v,t) = e-"*^^., (9) 



where the quark Dirac spinors have the form 



By using Eqs.(6)-(10), we obtain 



^=( ■^"H.'lx- (10) 



du(r) 
dr 



^{e + ga{r))v{r), (11) 



^ = --vir) + ie~gair))u{r), (12) 

dr r 

—Y + -—r^ --r^ Ng{u^{r) - v\r)). (13) 

dr^ r dr da 

The quark functions should satisfy the normalization condition 

4tt I r^{u^{r)+v^{r))dr = 1. (14) 

The number of quarks is A^ = 3 for baryons and N — 2 for mesons. In the following, our discussions are constrained 
in the case of TV = 3. These equations are subject to the boundary conditions which follow from the requirement of 
finite energy: 

.(0)^0,^^0, (15) 

dr 

u{oo) — 0, cr(oo) = a^. (16) 

If we consider N quarks in the lowest mode with energy e, the total energy of the system is 



E = Ne + An / r^ 



1 /dCT^ ^ 



2V*' ■"<-> 



dr. (17) 



The model has four adjustable parameters g, a, b, c which can be chosen to fit various baryon properites, such 
as masses, charge radii and magnetic moments. Let's set R, fip and gA/ 9v be the proton charge radius, the proton 
magnetic moment and the ratio of axial-vector to vector coupling respectively, then they are given by 

/•OO 

i?2 = 47r / r'^{u'^{r) + v^{r))dr, (18) 

/ip = — / r^u{r)v{r)dr, (19) 

•J Jo 



9a 207r [ 2t 2, \ ^2, N\ . 
— = — - / r [u [r)- -V {r))dr. 
gv i Ja 3 



(20) 



As in any simple quark model, the neutron magnetic moment is — | that of the proton, and its charge radius is zero. 
Once the solutions to the above equations are obtained, one can calculate these physical quantities pertaining to 
the three-quark system, which have been measured experimentally. In Ref. [10, 15, 16], a wide range of parameters 
have been used to calculate the quantities above. Hereafter we take one set of parameters a = 17.70fm~^, b = 
— 1457 Afm^^, c — 20000 and g = 12.16 to study the temperature and density dependence of the soliton solution. 
It has been proved in Ref. [10, 15, 16] that this set of parameters can describe the properties of nucleon at zero 
temperature successfully. 



III. THE ONE-LOOP EFFECTIVE POTENTIAL 

A convenient framework of studying phase transitions is the thermal field theory. Within this framework, the finite 
temperature effective potential is an important and useful theoretical tool. In this section, in order to investigate the 
temperature and the chemical potential dependence of the Friedberg-Lee model, we summarize the relevant results of 
Dolan and Jackiw[21], since the model we use is almost identical to the model studied by those authors. 

The finite temperature partition function Z{J; (3) is given by 

^u;a^ "-'-"°;[;ji;^-'"' . (21) 

where H is the Hamiltonian obtained from the Lagrangian(l) and (3 is the inverse of the temperature (/? = T~^). 
Texp(i J • a) denotes the time-ordered exponential of 

.,/..../...,;(.,.(.), ,22) 

where, in the imaginary-time formalism, x^ runs from to —ij3. The a field is required to satisfy periodic boundary 
conditions in time; the ip fields are antipcriodic in time. The generating function for the connected Green functions 
is given in terms of Z{J; (3) by 

W{J;l3)^-i\nZ{J:l3). (23) 

The effective action is obtained via a Legendre transformation of M^(J;/3), 

r(CT; /?) = VK( J; 13)- I d^xaJ{x), (24) 

where the average value of the a field is 

Here we study a uniform system in which a does not depend on the space-time coordinates. The generalized effective 
potential for such a system is defined by 

V{a;P)^-n-'r{a;f3), (26) 

where an infinite volume factor CI arises from space-time integrations and it is customary to introduce the generalized 
effective potential V(W;P). In the thermodynamic language, the thermal effective potential has the meaning of the 
free energy density and it is related to the thermodynamic pressure through the equation 

p^-V{a;f3). (27) 

The effective potential can be calculated to one-loop order by using the methods of Dolan and Jackiw[21]. However, 
it is pointed out in Ref[9, 22] that, as an approximation, all a quantum loop diagrams may be ignored due to the 
fact that a is only a phcnomenological field describing the long-range collective effects of QCD, and its short-wave 
components do not exist in reality. Therefore for the rest of this discussion we shall ignore quantum corrections and 
concentrate on those induced by finite temperature and density effects. 

In order to study the temperature and density dependence of the bag constant and the features of the deconfinement 
phase transitions, we first calculate the effective potential at finite temperature and density. Using the techniques 
of Dolan and Jackiw[21], the one- loop contribution to the effective potential of the Friedberg-Lee model at finite 
temperature and density is of the form 

V{a; 13, y) - U{a) + Vsio- /3) -f Vf{g- /3, m), (28) 

where U{(j) is the classical potential of the Lagrangian. Vb(ct;/3) is the finite temperature contributions from bo- 
son one-loop diagrams, and Vf{<j; P, ^) is the finite temperature and density contributions from fermion one-loop 
diagrams[19][21]. For simplicity, we have substituted (j{T) for (j{T). These contribute the following terms in the 
potential[21] 

1 POD ^^^ 

^B(a;/3) = 2^^^ / dxxHn[l - e-V^"^^^^) , (29) 
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FIG. 1: The one-loop Effective potential V{a; l3) as a function of a at T = lOMeV, T = 80MeV,T = 106.6MeV, T = 129MeV 
and T = 140MeV when the chemical potential ^ = 0. The critical temperatures is set at Tc — 106.6MeV when two minima 
are euqal. 



VFia;P,n) 



-12 






27r2/34 



dxx In I 1 



-{y/x^+l3^ml-0,i) 



(30) 



where the minus sign is the consequence of Fermi-Dirac statistics, mg 
field a and the quark field, respectively: 



and mq are the effective masses of the scalar 



a + ba{T) + ^a\T). 



(31) 
(32) 



We fix rn^ by taking its value at the physical vacuum state[15]. In this paper, the bag constant B(T,ii) is defined 
as the difference between the vacua of the effective potential inside and outside the solion bag. This means that for 
T < Tc (or /i < /ic)i B{T, ii) is the difference between the effective potential values at the perturbative vacuum state 
and values at the physical vacuum state 



BiT,^i)^V{ao;P,fl)~Via,■,P,^l). 



(33) 



For T > Tc (or /z > nc), the bag constant is zero due to the fact that the vacua inside and outside the soliton bag 
are equal. This will be analyzed in detail in next section. From Eqs.(28), (31) and (32), we can numerically solve the 
effective potential V{(7; (3, /i) for different temperatures and densities. 

When the chemical potential is zero, the effective potential at different temperatures is illustrated in Fig.l. It 
can be seen from Fig.l that there exist two particular temperatures. One is that the effective potential exhibits two 
degenerate minima at Tc — 106.6MeV which is defined as the critical temperature, the other is that the second 
minimum of the potential at cr ~ CTi, disappears at a higher temperature T ~ 129MeV. For low temperatures the 
absolute minimum of V{a; (3) lies close to (t„, and there is another minimum at Uq. The physical vacuum state at (j„ 
is stable, and correspondingly quarks are in confinement. As the temperature increases the second minimum of the 
potential at a^ decreases relative to the first one. At the critical temperature Tc — 106.6MeV, the potentials at the 
two minima are equal. The physical vacuum becomes unstable and soliton solutions are going to disappear. 

When fixing temperature T = 50MeV, in Fig. 2 we plot the one-loop Effective potential y(i7;/3,/x) as a function 
of a at difference chemical potentials ^J. = lOOMeV, n = 223.1MeV, fi = 300MeV and /i = 354MeV . From Fig.2 , 
there are two specific chemical potentials, one corresponds to the critical chemical potential iic — 223.1MeV where 
the effective potentials at the two vacuum states are equal, the other is a higher chemical potential /i ~ 354MeV 
at which the chemical potential increases and the second minimum of the potential at ct ~ cr^, disappears. For low 
chemical potentials the absolute minimum oi V{a; /3, n) lies close to ay, and there is another relative minimum at aQ. 
The physical vacuum state at ay is stable and correspondingly quarks are in confinement. As the chemical potential 
increases the second minimum of the potential at CTo decreases relative to the first one. At the critical chemical 
potential fic — 223.1MeV, the potentials at the two minima are equal. When the chemical potential increases further, 
the physical vacuum becomes unstable and the soliton solutions tend to melt away. 
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FIG. 2: The temperature and density dependence of the one-loop effective potential V(a\ (3, p) when fixing temperature T at 
50 MeV. At low chemical potential, there are two minima and one of which is going to disappear when the chemical potential 
reaching a high value ji — 354MeV. The critical chemical potential is set at fic — 223. f MeV when two minima are euqal. 
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FIG. 3: The bag constant B{T) as functions of T when the chemical potential is zero. 
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FIG. 4: The bag constant B(T, /i) as functions of /i when the temperature is fixed at T = 50MeV. 
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FIG. 5: The phase diagrams for the Friedberg-Lee model in the (T — /i) plane. Along this line the effective potential at the 
two vacuum states are equal, the physical vacuum becomes unstable and the soliton solutions tend to melt away. 

From the above effective potential, it is convenient to investigate the temperature and the chemical potential 
dependence of the bag constant. For the case of finite temperature and zero density, the bag constant B{T) in the the 
Friedberg-Lee model is defined as the difference between the vacua of the effective potential inside and outside the 
solion bag. This means that for T < Tc, the bag constant is defined as the difference between the effective potential 
at the perturbative vacuum state and the effective potential at the physical vacuum state. 



BiT)^V{ao:l3)-Via,;P). 



(34) 



For T > Tc, B{T) — 0. From Eq.(34) we illustrate the temperature dependence of the bag constant B{T) in Fig3. 
Since B{T) — when T is lager than Tq, there is no bag constant to provide the dynamical mechanism to confine 
the quarks in the bag. Moreover, there is no soliton-like solution in the model when the temperature is above Tc 
due to the effective potential and vacuum structure, there only exists the damping oscillation solution which can not 
produce a mechanism to confine the quarks in a small region, we will investigate the damping oscillation next section 
in details. Therefore, as the temperature is above Tc, the confinement of the quarks is removed completely. Similarly, 
for the case of fixing temperature at T = 50MeV, the bag constant at different densities is defined as 



B(T,^i) = v{<Jo■,p,^i)-v(a,■,|3,^i). 



(35) 



And for fi > fxc, B{T,ij) = 0. The result is shown in Fig. 4. In this case, we can see the deconfinement phase 
transitions begin to take place at the chemical potential /ic- 

In the end of this section, by using the similar discussion in Ref.[23], the phase diagram in the T — fi plane 
calculated for the Friedberg-Lee model is shown in Fig. 5. The critical line corresponds to the phase transition line 
between confined and deconfincd phases. At the critical T and /x the two vacumm states, namely the perturbative 
vacumm state and the physical vacumm state, appear in the same energy. This is treated as the sign of the first order 
phase transition[15][17] [23]. Our phase diagram shown in Fig. 5, based on the phenomenological Friedberg-Lee model, 
predicts a first order deconfinement phase transition for the full phase diagram since there are two equal minima in the 
effective potential separated by the barrier. This result is different from the predictions based on lattice gauge theory 
where more complicated phase diagram on the T and /i plane has been obtained, giving more fruitful behaviors of the 
QCD phase [24] [25]. For example, at finite temperature and zero chemical potential, the phase transition could be of 
a second-order, moreover, it is predicted that there may exist a critical point along the critical line, at lower baryon 
density there might be a rapid crossover, while at higher baryon density, the first order phase transition is obtained. 
Finally, we would like to point out that the assumed phenomenological potential in the Friedberg-Lee model should 
be the reason for these differences in the behavior of phase diagram. 

IV. THE SOLITON SOLUTIONS AT FINITE TEMPERATURE AND DENSITY 



Field models such as those described above can be modified so as to allow for the effect of a thermal background, 
with temperature T, just by replacing the relevant classical potential function U{cr) by an appropriately modified 
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FIG. 6: The solutions for the set of coupled nonlinear differential equations (11), (12) and (36) for different temperatures when 
T <Tc and n = OMeV. 
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FIG. 7: The solutions for the set of coupled nonlinear differential equations (11), (12) and (36) for different temperatures when 
T > To and ^ = OMeV. 



temperature-dependent effective potential V(cr; /3) [19] [26] [27]. In this work, in order to study the effects of the tem- 
perature and density on the soHton solutions of the Fricdberg-Lcc model, the classical potential function U{a) should 
be replaced by the temperature and density dependent effective potential V{a; (3, fx). 

As in the zero temperature and zero chemical potential case, for finite temperature and density soliton solutions will 
take the same form as the one at zero temperature and zero chemical potential. But the a field should be determined 
by the equation of motion: 



dV(r) 2dcr(r) dV{a;l3,^i) 
dr"^ r dr da 



Ng{u^{r) -w^(r)). 



(36) 



where the temperature and density dependent effective potential V{a;(3,fi) is defined in Eq.(28). And the quark 
functions should also satisfy the normalization condition 



47r / r^{u^{r) + v {r))dr ^ 1 



(37) 



The functions cr(r), u{r) and v{r) also satisfy the boundary conditions following from the requirement of finite 
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energy: 



v(0) 
u{oo) 



rfcr(O) 
dr 



(38) 
(39) 



However, the situation is changed according to (T(r) as r — > oo. When T < Tc (or fi < ^c)i in order to satisfy the 
requirement of finite energy of the sohtion (or other topological defects), as r ^ oo, a{r) should be equal to a^, where 
the potential V{a) has an absolute minimum. For T > Tq (or /i > /ic)i the physical vacuum becomes unstable, and 
the stable vacuum is the perturbative vacuum which is the absolute minimum of the effective potential. Therefore, 
to satisfy the requirement of finite energy of the solition at T > Tc (or fi > lie), we should take the asymptotic value 
(vacuum values) cr —^ as r ^ oo. Based on above analysis, one obtains the following boundary condition for the 
function u(r) as r ^ oo: 



a = ao, for 



for T <Tc 
Tc <T or 



or ^ < iJc, 
He < M- 



(40) 
(41) 



In order to investigate the behaviors of soliton solutions at finite temperature and density, the set of coupled nonlinear 
differential equations (11), (12) and (36) should be solved numerically with the boundary conditions Eqs.(38),(39),(40) 
and (41). In the following discussion, we will focus on two cases, one is that the temperature is T < Tc as the chemical 
potential is zero, the other is that the chemical potential is /i < iic by fixing the temperature at T = 50MeV. 

For the first case, in Fig. 6, we plot the soliton solutions by taking the temperatures as T = OMeV, 60MeV, QOMeV 
and 106.6MeV. From Fig. 6, we can see that with the temperature increasing, (Jt,(/3) is nearly constant near ay, 
while (7(0) is changed dramatically, and there exist stable soliton bags for T < Tc- Also it can be seen in Fig. 6 the 
radius of a stable soliton bag increases as temperature increases. When the temperature is Tc < T, in Fig. 7, we 
plot the solutions by taking the temperatures as T = 11 OMeV, 120MeV, 123MeV and 123MeV. As the temperature 
is higher than the critical temperature Tc, the bag constant is zero, the physical vacuum becomes unstable, and the 
perturbative vacuum state is stable, the boundary condition for the function a is changed accordingly, then the shapes 
of solutions are very different from that of T < Tc- With the increase of temperature, the soliton solutions tend to 
disappear. From Fig. 7 we can see that, unlike the conventional soliton solutions, here the solutions have the damping 
oscillation, and we can not find soliton solution anymore. 

Following the discussions of Ref. [19], the radius of a stable soliton bag increases as temperature increases. From 
Eq.(18), in Fig. 8 we plot the variation of radius as functions of temperature, which reveals that the radius of the 
soliton bag does increase with increasing temperature. Once the temperature is larger than the critical temperature 
Tc, the physical vacuum becomes unstable while the perturbative vacuum state is stable. Moreover, just as above 
mentioned, there is no soliton solution when the temperature crosses over the critical temperature, this means that 
the stable soliton bag is going to melt away and disappear. Therefore there is no definite radius for T > Tc- 

Similar discussion can be applied to the second case, when the temperature is fixed, the soliton solutions at different 
chemical potentials are illustrated in Fig. 9, where the chemical potentials are taken as /i = lOOMeV, 160MeV and 
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FIG. 9: The solutions for the set of coupled nonlinear differential equations (11), (12) and (36) for different chemical potential 
when n< He and T = 50MeV. 
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FIG. 10: Radius of the soliton bag as functions of chemical potential when /i < /ic and T = 50McV. R02 is the radius of the 
soliton bag for T = 50MeV and h = OMeV. 



223.1McV for a fixed temperature T — 50MeV. With the chemical potential increasing, the radius of a stable soliton 
bag increases too, we plot the radius of the sohton bag versus chemical potential at a fixing temperature in Fig. 10. 
We can see clearly that the radius of the soliton bag increases continuously with the increase of chemical potential. 
When the chemical potential is larger than the critical chemical potential no, due to the fact that the physical vacuum 
becomes unstable and the perturbative vacuum state is stable, there is no definite radius and the bag does not exist 
anymore. In order to show that there is no stable soliton bag when the chemical potential is above the critical chemical 
potential, the set of coupled nonlinear differential equations (11), (12) and (36) should be solved numerically with 
the boundary conditions Eqs.(38),(39) and (41) for y^c < M by fixing the temperature at T = 50MeV. We plot the 
solutions when the chemical potential is larger than the critical chemical potential by taking the chemical potential 
as 11 = 225MeV, 240MeV, 260MeV, 270MeV and 277MeV in Fig.ll. When the chemical potential is higher than the 
critical chemical potential fxc, the bag constant is zero, the physical vacuum becomes unstable, and the perturbative 
vacuum state is stable, accordingly the boundary condition for the function a is changed, then the shapes of solutions 
are very different from that of /i < iic- With the increase of chemical potential, the soliton solutions tend to disappear. 
From Fig.ll it can be seen that the solutions have the damping oscillation and there is no soliton solution. 

As mentioned in Sect. Ill, B{T) is defined as the difference between the vacua inside and outside the soliton bag. For 
T <Tc and /u = OMeV, the vacuum inside the soliton bag is the metastable vacuum, the vacuum outside the soliton 
bag is the real physical vacuum. This can be seen in Fig. 6. As the temperature is above the critical temperature 
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FIG. 11: The solutions for the set of coupled nonlinear differential equations (11), (12) and (36) for different chemical potential 
when // > //c by fixing the temperature T = 50MeV. 

Tc, the metastable vacuum becomes the absolute one. The behavior of cr(r) changes dramatically when T > Tc^ i.e. 
a(r) is very close to zero at any r, this can be seen from Fig. 7. Therefore, the values of the effective potential are 
always close to their stable minimum at a = 0. The state at cr ~ '^v{P) can no longer be realized in this case. This 
means that the bag constant B — 0. Based on above analysis, there is no soliton solution and the bag constant B{T) 
is zero for T > Tc, so there exists no more mechanism in Friedberg-Lee model to confine the quarks, and the quarks 
are to be deconfined when T > Tc- Accordingly, for /i < ^c and T = 50MeV, the vacuum inside the soliton bag 
is the metastable vacuum, the vacuum outside the soliton bag is the real physical vacuum shown in Fig. 9. As the 
chemical potential is above the critical chemical potential /ip, the metastable vacuum becomes the absolute one. It 
can be seen from Fig. 11, the behavior of a{r) changes dramatically when /x > /dc, i.e. cr(r) is very close to zero at 
any r. Therefore, the values of the effective potential are always close to its stable minimum at cr = 0. The state 
at CT ~ (Tt,(/3) can not exist in this case, and the bag constant B — 0, hence there is no soliton solution and the bag 
constant B{T) is zero for ^ > /ic, no more mechanism in soliton bag model can be applied to confine the quarks, 
then the quarks are to be decofined in high chemical potential. 



V. SUMMARY AND DISCUSSION 



In this paper, we have studied the deconfinement phase transtion of the Friedberg-Lee model at finite temperature 
and chemical potential. Our studies are mainly concentrated on the effective potential of the Friedberg-Lee model 
and the soliton solutions at different temperatures and chemical potentials. We have gotten the bag constant functions 
B{T) and B{T, fx) for finite temperature and density respectively. It is shown that the two minima of the potential at 
zero temperature and chemical potential will be equal at a certain temperature T ~ 106.6MeV or a certain chemical 
potential n ~ 223.1MeV, which are to be defined as the critical temperature Tc and the critical chemical potential 
^c accordingly. 

When the temperature T > Tc, the original perturbative vacuum state becomes stable and the original physical 
vacuum state becomes a metastable one. Because cr(r) is very close to zero at any r for T > Tc, the effective potential 
always takes its value at the absolute minimum V{<t ~ 0). This gives B{T) = for T > Tc- As the temperature 
approaches another higher temperature T ~ 129MeV, the effective potential has a unique minimum. Similarly, for 
/i > fic, the original perturbative vacuum state becomes stable and the original physical vacuum state becomes a 
metastable one. Due to the fact that cr(r) is very close to zero at any r for /i > fic, the effective potential always 
takes its value at the absolute minimum V{(7 ^ 0). This makes B{T,ii) — for fi > iic- As the chemical potential 
approaches another higher chemical potential /i ~ 354MeV, the effective potential has a unique minimum. Our results 
are qualitatively similar to that obtained in conventional soliton bag model at finite temperature[15, 16, 17, 19]. 

In Friedberg-Lee model, the confinement of the quark requires the existence of the soliton solution, and the soliton 
solution depends on the effective potential at finite temperature and density. In order to investigate the behavior of 
the soliton solution at finite temperature and chemical potential, we numerically solve the set of coupled nonlinear 
differential equations. Our results show that when T < Tc (or ji < ^c), there exist the stable soliton solutions in 
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the model, however, when the temperature is higher than the critical temperature Tq — 106.6MeV (or the critical 
chemical potential jiq — 223.1MeV), there is only damping oscillation solutions and no soliton-like solution exists, 
and the quark can not be confined by such solutions anymore, then the confinement of quarks are removed and the 
deconfined phase transition takes place. We also obtain that the radius of the soliton bag does increase when the 
temperature or density increases. At T = Tq (or /x = fic) the soliton bag disappears. 

Since the soliton bag model has four adjustable parameters, the numerical results of the critical temperature and 
the critical chemical potential of the deconfinement phase are parameter-dependent. Moreover in this work, we have 
only studied the single nucleon properties at finite temperature and chemical potential and have not considered the 
interactions of the nucleons. In order to obtain a reasonable description of nuclear interactions [29] [30] [31], we should 
consider the modifications of the Friedberg-Lee model to include explicit meson degrees of freedom coupled linearly 
to the quarks, which have been studied recently in Ref. [32] at zero temperature, it is of interest to extend their work 
to finite temperature and density, and discuss the soliton solutions at different temperatures and chemical potentials. 
All these works are in progress. 
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